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We study energy-transport efficiency in light-harvesting planar and 3D complexes of two-level 
atomic quantum systems, embedded in a common thermal blackbody radiation. We show that the 
collective non-local dissipation induced by the thermal bath plays a fundamental role in energy 
transport. It gives rise to a dramatic enhancement of the energy-transport efficiency, which may 
largely overcome 100%. This effect, which improves the understanding of transport phenomena 
in experimentally relevant complexes, suggests a particularly promising mechanism for quantum 
energy management. 


INTRODUCTION 

Transport of energy is a phenomenon of the utmost im¬ 
portance in several domains of physics. Recent advance¬ 
ments in manipulations of systems approaching quantum 
scales have further triggered research on transport be¬ 
tween single quantum entities m- This interest has a 
two-fold reason: first, the understanding of energy fluxes 
is necessary for controlling micro- and nanoscale devices 
[B H 13; second, purely quantum effects allow for en¬ 
ergy management not achievable in classical scenarios Et 
[8]. Probably the most celebrated of these effects is the 
role of quantum coherence in photosynthetic efficiency 
0110]. Stimulated by these results, quantum excita¬ 
tion transport in networks of two-level systems (TLSs) 
has recently received an enormous deal of attention m 
In this context, the interaction of TLSs with vari¬ 
ous types of environments has been largely investigated 
[mils US- In particular, the effect of thermal black- 
body radiation has been studied using phenomenological 
models only including local noise acting on each atom 
P^SHTT] . which typically deteriorates transport as tem¬ 
perature increases. This approach, neglecting non-local 
effects induced by the radiation, led the attention mostly 
on low-temperature transport, where non-locality is in¬ 
deed marginal. 

In this Letter, through a microscopic model, we ex¬ 
plore excitation transport along TLSs in a common ther¬ 
mal blackbody radiation at arbitrary temperature, fully 
taking into account radiation-induced non-local effects. 
We show that the collective non-local dissipation due to 
the correlations of the thermal radiation may have a dra¬ 
matic effect on transport: unexpectedly, transport effi¬ 
ciency may increase with temperature up to values far 
beyond 100%. 


PHYSICAL SYSTEM 

We consider N identical two-level quantum emitters 
(hereby referred to as atoms) of frequency uJa with ground 


(excited) state jO^) (|li)), embedded in a blackbody radi¬ 
ation at temperature Tb- They couple to the electromag¬ 
netic field through their dipole moment operators, having 
ground-excited matrix elements all equal to d. Without 
external perturbations, the atomic system equilibrates at 
Tb- Suppose now that two of these atoms are externally 
perturbed by pumping excitations into one and extract¬ 
ing energy from another (see sketch in Fig. [^. 

In the weak-coupling limit, the dynamics of the atomic 
density matrix p is described by the Markovian master 
equation (ME) [18] 

P — ^ \_Htot, p] + F)\qcP -|- DyiIP D[yiP T -DoutP: (1) 

where fftot = Ha + H^, Ha = EiHi = 

hu>a Ei F is the atomic free Hamiltonian and hF = 
hKijCjFcF^ the field-induced atomic dipolar coupling, 

with ^ the raising (lowering) operator of atom i. 

D\ocP = E (Npcrr ~\Fi F > p}) 

1 1 

+ (1 + n)7o (cr^ paf - -{af , />}) | 

describes standard thermal dissipation and 70 = 

d‘^uj^/{3h7reoc^). The photon number n = n{TB,00a) = 
l/[ex.p{huJa/kBTB) — 1 ], encompasses the fundamental 
dependence of atomic dynamics on T^. Moreover 

Dn\p = E Fpp = EI {F pf ~ F > f^}) 

i<j i<j ^ 

+ (1 + n)jij (a~paj - ^ {(7+, p}) + h.c. | (3) 

describes two-body non-local dissipation, characterising 
energy exchanges between the field and any two-atom 
subsystems. These terms, neglected in phenomenologi¬ 
cal descriptions of transport based on atom-field coupling 
naiiaiin], are a natural consequence of non-local field 
self-correlations in a microscopic derivation of the ME 
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[iHiiin]- Their effects are a manifestation of the quan¬ 
tumness of atom-field coupling: they have indeed been 
studied with respect to bipartite entanglement of spins 
in a common environment [2TI423] . The rates involved in 
the ME 0 are analytical functions of interatomic dis¬ 
tances, atomic dipoles and atomic frequency m- In the 
case of identical and real dipoles d considered here, they 
result 


Incoherent pumping on the atom labeled as p is ac¬ 
counted for by the term 

DinP = 7in ^cr+ pa- - i cr+ , (8) 

Another external operator (the sink) extracts energy lo¬ 
cally from the atom labeled as e, through the process 


=70 E (W") (4) 

n=l,2,3 


-^outP — Tout { 


1 

2 





(9) 


where d = d/|d|, [d]i is the component along the line 
joining the two atoms and [d]( 2 , 3 ) the components along 
directions perpendicular to each other and to the line 
joining the atoms. Being = Vj — the spatial sepa- 



FIG. 1. Stationary transport efficiency Xstat versus Tg for 
three planar configurations of 5 atoms on a ring. Atoms are 
either placed on the vertices of a regular pentagon (RP, blue 
points and dot-dashed line), one of them is displaced by an 
angle O 2 = —0.37 rad (clockwise-hopping CH, green points 
and dashed line), or two of them are displaced in a symmetric 
way with Op = —O 2 = 0.37rad (no hopping NH, orange points 
and solid line). 


ration of two atoms, Vij = Vij = y^, coefficients 

(n) 

alj are 

“T = ^ ( ) ’ (5) 

“if = “if = ^ (^ii cosfij + {fl - 1) sinfy). (6) 
Finally, the dipole-dipole interaction strength is 


In general, a Markovian dissipative term D produces 
an energy flux Tr (iT^Dp), Hs being the free Hamiltonian 
of the open system [24]. According to this definition, 
energy flowing out of an open system is negative. The 
rates of energy pumping P{t) and extraction E{t) are 
given by 

P{t) = An/9(i)) , (10) 

E{t) = ^^UE{t) = -Tv{HaDoutp{t)), (11) 

Up(E){t) being the total energy pumped in (extracted 
from) the atomic system up to time t. By construction 
P{t)^E{t) > 0, and depend on Tp through the atomic 
state p{t). Due to thermal excitations, energy is ex¬ 
tracted from e also in the absence of pumping (yin = 0): 
the sink, through atom e, absorbs thermal photons from 
the field at a rate Eq. To evaluate how effective the 
transport of excitations from p to e is, one can analyse 
the contribution to extracted flux due to pumping only, 
as expressed by the transport efficiency 


x{t) 


E{t)-Eo{t) 

P{t) 


( 12 ) 


When x(t) > 1, E > Eq P, i.e., the extracted ex¬ 
citations are more than the sum of the thermal and the 
pumped ones: as will be shown in the following, efficiency 
will largely overcome 100% thanks to the combination of 
thermal noise and non-local effects. The definition of y 
does not include energy fluxes between atoms and the 
blackbody radiation and, as such, y is not upper-limited 
by 1. It is therefore not a thermodynamic efficiency but 
rather characterises how efficiently pumping affects ex¬ 
traction, as these are the only two quantities on which 
an external control is possible. 

At stationarity, all the quantities involved in Eq. (12) 
become time-independent, hence both Up and Up di¬ 
verge linearly in t and 


Aij = -^70 [2(d ■ rij)^f{fij) + (1 - (d ■ , 

. COS X-h X sin X , . (x^ — 1) cosx — xsinx 

/(^) = -^3-’ = -^3-• 

( 7 ) 


Xstat 


r/'stat r/'stat 

E _ Eq 

[/stat 


(13) 


is the fraction of total pumped energy one recovers at the 
extraction site. 
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Following Eq. 0 ’ atoms have two ways of affecting 
each other’s internal energy: coherent interaction Hi 
and collective dissipative exchanges Dfj with the envi¬ 
ronment. The associated energy fluxes are = 

p{t)]) (flux of energy from i to j) and 

Qfj{t) = Ti(^Hjp{t)) (energy flowing non-locally 
from the field to each of the atoms i,jf, i.e., two pho¬ 
tons are emitted or absorbed simultaneously by the two 
atoms). Explicit evaluation gives 

Qfj (t) = - {crt crj + cr ■) !■), (14) 

= ihWaAij{erf (j- - (cr+ aj )t). (15) 

The correlators of emission and absorption processes 
satisfy {crterj + (cr+cr“)’*') = 2 Re(c*-^’(t)) and {crterj - 
{(jf (jJ^) = 2 ilm(c*-^(t)), {t) being the quantum co¬ 

herence between the states \0ilj) and |liOj) in pij = 
Tia^ijp. We remark that dissipation plays here an ac¬ 
tive role in transporting excitations coherently [Eq. 0 ], 
contrarily to previous models where its effect mostly con¬ 
sisted in localising excitations by suppressing coherences. 

TRANSPORT IN PLANAR CONFIGURATIONS 

We apply now this theory to transport through five- 
atom rings of radius r = 0.4//m, as in Fig. Pla¬ 
nar ring configurations have been studied as a simpler 
prototype of light-harvesting physics El- The atomic 
frequency is Ua = 10 ^^ rad/s {huJa — b .6 • 10“^ eV) and 
dipoles lie along the axis perpendicular to the ring. These 
parameters are in accordance with some typical values 
found in, e.g., quantum dots m, hyperfine structure of 
real atoms and molecular vibrational states. Moreover, 
although not usually met in biological light-harvesting 
complexes [26], these parameters will permit us to high¬ 
light more clearly the general importance of the non-local 
dissipation in transport phenomena. The pumping rate 
is 7 in = 10 “^ 7 o, to simulate low-intensity light impinging 
on the system, while the extraction rate 7 out = 10 ^ 70? 
in accordance with typical sink rates in light harvesting 
complexes [26l[27]. Thanks to this choice, all the term in 
the ME are proportional to 70 . In what follows we will 
therefore work in units of 70 . Pumping and extraction are 
exerted as shown in Eig. All the results reported here¬ 
after are obtained using the open-source QuTiP package 

m- 

Three different geometrical dispositions are considered. 
Eirst, atoms are placed on the vertices of a regular pen¬ 
tagon (RP), as shown in Eig. (blue points). In this 
case two paths, clock- and counterclockwise, exist for ex¬ 
citations to travel to atom e. If is high enough, it 
is known that field-induced effects can create multipar¬ 
tite correlations in symmetric atomic configurations [29| : 
excitations get delocalised over the whole atomic system 


and transport efficiency is strongly reduced. Eig. shows 
that transport is indeed rapidly killed by thermal noise. 

The second configuration is obtained through a dis¬ 
placement of atom 2 from its regular position of an an¬ 
gle O 2 = —0.37rad (green dots in Eig. [^. We refer to 
this configuration as clockwise hopping (CH): the fast 
decay of Aij with the interatomic distance Vij induces a 
bias between the two hopping directions. Excitations are 
mostly transmitted along the path p ^ 5 ^ e with fewer 
destructive interference: Xstat for this configuration is de¬ 
picted in Eig. [^(dashed green line). In accordance with 
hopping physics, local noise disturbs excitation transmis¬ 
sion, continuously reducing its efficiency with increasing 
Tb. 

Einally, also the second pathway is killed by displac¬ 
ing the atom p closer to atom 2, such that a symmet¬ 
ric configuration is obtained (orange points in Eig. [^. 
A natural expectation is that transport becomes hardly 
possible due to the vanishing hopping amplitude (no¬ 
hopping configuration NH). Nonetheless, Xstat shows a 
pronounced maximum in Tg, at efficiencies higher than 
300%: remarkably when hopping is killed, the increase 
in extracted excitations becomes larger than pumping. 
There being no hopping, such an enormous gain of exci¬ 
tations origins from by which atoms can affect each 
other’s internal energy through collective field-induced 
dissipation. 

To better understand this mechanism, the atomic sys¬ 
tem can be divided into three subsets: the pair {p, 2} 
where excitations arrive, the pair {3,5} which trans¬ 
mits these excitations and the monomer {e} in contact 
with the sink, as shown in the inset of Eig. m Only 
a small fraction of excitations jumps directly from p 
to e, the largest part travelling along the main path 
{p, 2} ^ {3,5} ^ {e}. Since local heat fluxes can not 



Tb(K) 


FIG. 2. xSat (orange solid line), xSat (green dashed line) 
and the principal contributions to all stages of the main path. 
The first stage {p, 2} ^ {3, 5} in CH and NH configurations 
happens through different physical effects: hopping 77 for CH, 
non-local dissipation ly for NH. 
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affect transport due to the Markovian structure of the en¬ 
vironment, the only possible mechanisms responsible for 
excitation transmission between atoms are hopping and 
non-local dissipation. To analyse them, we introduce 



a=a,b 

K,=i,j 



a=a,b 

K,=i,j 


representing the fraction of pumped excitations travelling 
to the subset {a, b} through either hopping {r]) from or 
non-local dissipation (z/) with each atom of the subset 
{i, jf}. Fig. [^shows, for both NH and CH configurations, 
the quantities Xstat and versus Tb- It is clear that 
the mechanism responsible for the last stage {3, 5} ^ {e} 
of the main path is hopping: the curves (solid 

orange and dotted green lines) and (dotted black 

and dashed blue lines) are remarkably similar, showing 
that almost all the energy extracted from e arrives there 
through hopping from atoms 3 and 5. 

A striking difference appears however in the first stage 
{p, 2 } ^ {3, 5}: in CH case, it is again mostly due to hop¬ 
ping from pair {p, 2} to pair {3,5}, as shown by ^^ 2^35 
(brown dot-dot-dashed line). Non-local transport ^^^35 
(not shown here) is indeed negligible. On the other hand, 
the stage {p, 2} ^ {3, 5} in NH setup has a large con¬ 
tribution from collective non-local dissipation as high¬ 
lighted by the non-local energy exchanges ^^^ 2^35 (red dot- 
dashed line): the pumping of excitations on p triggers the 
electromagnetic field to inject photons into both atomic 
pairs. Note that, as a consequence, excitations received 
by {3, 5} can exceed the ones pumped into p, the differ¬ 
ence being supplied by the field: ^^^ 2^35 accounts almost 
exactly for the additional parts of excitations travelling 
to e. Energy arriving on {3, 5} completes then its path 
by hopping to e. Following Eq. (14), non-local energy 
flux is proportional to the real part of the quantum co¬ 
herence between two atoms. Transport in NH configu¬ 
ration is thus a dissipation-induced coherent effect. The 
importance of the symmetric positions of atoms p and 2 
in creating a pair is proved by the dependence of Xstat 
on the angular displacement Op of p, reported in Fig. 
There, Op = 0 corresponds to CH configuration. Note 
that the role of symmetry on quantum correlations and 
quantum transport was already pointed out in different 
contexts [29ll3^ . 

The reason why these non-local effects do not con¬ 
tribute when hopping is significant is twofold: first, when 
hopping is active there can be no isolated subset of atoms 
where pumping is exerted, while such subset is present 
in NH configuration (pair {p, 2}). This subset, having a 
much larger excited-state population than the rest of the 
system due to pumping on it, drives the collective atomic 
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Tb(K) 

FIG. 3. Xstat versus Tb and angular displacement Op of p, 
with O 2 = —0.37 rad. The dashed white line is the contour 
Xstat = 1- The horizontal dot-dashed line represents the NH 
configuration (Op = —O 2 = 0.37rad), where transport effi¬ 
ciency is maximised. 



state further away from a purely thermal state, thus al¬ 
lowing the establishing of steady quantum coherence: we 
develop in the next section a simple toy model able to 
grasp this mechanism. The second reason comes from 
the timescales of different dynamical processes: Fig. 
shows the time- and Tg-dependence of x(E [panel 

(a)], CH [panel (b)] and NH [panel (c)] configurations 
with initial state IV^o) = ( 8 ) 2=1 | 0 «)* 

Compare panels (b) and (c): transport through hop¬ 
ping (b) happens on much shorter timescales than the 
non-local one, as expectable since A^j ^ As long 
as hopping is active, non-locality cannot play any sub¬ 
stantial role: by the time non-locality becomes effective, 
excitations have already travelled through hopping. Ex¬ 
citations reach e in jot < 1 in CH setup and stationarity 
is attained for 70 1 < 10. On the other hand, be- 

comes non-negligible in yot > 10 , while stationarity is 
reached for yot > 600. Nonetheless, at suitable tempera¬ 
tures, yot 100 is enough for transport efficiency in NH 
configurations to go beyond 1 (white dashed line in panel 
(c)), fully entering the non-local amplification regime. 


THERMALLY-INDUCED COHERENCE 
AMPLIFICATION: A TOY MODEL 

The pairs {p, 2} and {3, 5} in NH configurations are 
almost not connected by hopping. Excitations pumped 
in {p, 2 } can thus hardly leave the pair, and the excited 
state populations of atoms p and 2 grow remarkably with 
respect to those in {3,5}. The effect is thus to create a 
gap in populations between two branches of the atomic 
system. This means that the system itself is no longer in 
the vicinity of thermal equilibrium and can thus sustain 
quantum coherence between its subparts. To understand 
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FIG. 4. Density plot of x(t) for (a): RP, (b): CH, (c): NH, versus time and Tb- The white dashed line in (a) represents the 
contour x{t) — 0? while in (c) it represents the contour x(t) = 1. Neither of the two contours appear in CH setup. Dynamics 
starts from li/^o) = (8)i=i |0*)- 


this on a more quantitative basis, consider a toy model 
of two atoms, having the same ME as Eq. 0 but with 
^ij = Tin = Tout = 0- Moreover, a population gap is 
imposed by two different blackbody reservoirs, one for 
each atom, such that n in Eqs. 0-0 is different for 
the two atoms. In particular, one atom is connected to 
a hot reservoir with n = Uh (simulating the effect of 
pumping), while the other one is connected to a reservoir 
with n = ub < Tih- Non-local effects are also assumed to 
happen at The atoms coherence c in the decoupled 
basis is always real and depends on and as 

/ X + n-h + nB){nB - Uh) 

where 

F{nB, Uh) = 7o (1 -h 2nB)(l -h 2n/i)(l -h ns -h nuf 

- 7^1 [l -I- ns(7 -I- 13nB -l- 8n^) -|- (19) 

-I- 4nftnB(3 -h 6nB -h 4n|) - n\]. 

Under the natural condition 70 > TnU \c{nB^Tih)\ shows 
a maximum as a function of at Despite its 

simplicity, then, this toy model gives already important 
hints to better understand the physics of non-local trans¬ 
port: driving a collective atomic state far from its Gibbs 
form by imposing a population gap between two parts 
generates real quantum coherence between them, with a 
sharp peak when an optimal gap is achieved. This, in 
turn, means a sharp peak in the non-local energy flux in¬ 
volving the two subparts. The symmetric displacement 
characterising NH configurations produces exactly such 
an optimal value of population gap. 


TRANSPORT OF EXCITATIONS IN 
FMO-INSPIRED CONFIGURATIONS 

We consider next a 3D setup inspired by the cele¬ 
brated Fenna-Matthews-Olson (EMO) complex [26], a 
pigment-protein complex used for photosynthesis. Light 
is collected and transmitted by 7 molecular complexes 
(schematically seen as TLSs) to a sink, where it is used for 
biological functions. Nowadays, its high light-transport 
efficiency is believed to have quantum origins |9j[26]. Its 
transport properties, as a function of temperature of a 
photonic bath, have been studied in [16], but non-local 
dissipation has there been neglected. 

The spatial configuration of the 7-site system we study 
is shown in Eig. [^ obtained from the EMO site configura¬ 
tion (given for instance in [31]) by rescaling each site-site 
distance by a factor 700. Typical lifetimes of excitations 
in EMO complexes are ^ 1 ns m. while extraction time 
can be of the order of 0.25 ps: we set Tout/To = 2.5 • 10^. 
The ratio Tin/To = 10“^ is compatible with the low- 
intensity light impinging on the molecule, having verified 
this value to be low enough for the qin-dependence of 
Xstat lo be flat. There being no preferential direction is 
this system, we choose the dipoles to lie along the x-axis 
of Eig.[g 

It is important to stress that, by this configuration, we 
do not intend to analyse the exact physics of transport in 
EMO complexes where transition frequencies, dipole ori¬ 
entations, intersite distances and environment are much 
different from the ones considered here [H [261 [27l|3l]. 
Nevertheless, since the relative weight of non-local effects 
over hopping is ultimately due to configurations of sites in 
a network, uncovering their role in a geometry analogous 
to the EMO might shed light on equivalent effect dur¬ 
ing photosynthesis and suggest how to enhance transport 
performances in artificial realisations. The EMO geome- 
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FIG. 5. Xstat versus Tb in FMO geometry, with site-site dis¬ 
tances scaled with respect to the FMO complex by a factor 
0.7 • 10^. The black dashed line represents the case = 0, 
the solid red line is the efficiency obtained when ^ij ^ 0. 

try shows indeed two sites (p and 2 in Fig. [^, relatively 
closer to each other than to the rest of the complex. 

Plots in Fig. show Xstat versus Tg in the absence 
(dashed black line) and in the presence (solid red line) 
of ^ij: the presence of non-local effects produces two in¬ 
teresting features: first, Xstat is always higher if ^ 0. 
Second, it introduces a non-monotonic dependence on Tg 
very similar to the one shown in Fig. suggesting an 
important contribution of non-local energy transfer. 

For completeness, we have also simulated the exact 
FMO complex, where absorption happens mostly around 
A = 800 nm (cj^ = 2.35 • 10^^ rad/s) [^, using the exact 
intersite distances and dipole orientations [31]. For this 
high frequency, n{TB,00a) does not significantly change 
around room temperature and the Tg-dependence of 
Xstat is flat. Nonetheless, non-local effects still play a sig¬ 
nificant role: Xstat - 0.8 when jij = 0, while Xstat - 0.88 
if jij 7 ^ 0, increasing transport efficiency of ^ 10%. 

CONCLUSIONS 

We showed that the energy transport in atomic quan¬ 
tum systems can be considerably affected by embedding 
them in an electromagnetic thermal bath. Non-local ef¬ 
fects induced by the thermal field, neglected in previous 
investigations, are shown to play an active role in coher¬ 
ent excitation transport and to boost its efficiency, which 
grows with temperature, up to values largely overcoming 
100% in planar configurations. This mechanism origi¬ 
nates from non-local self-correlations of electromagnetic 
radiation, which induce collective dissipation in atoms: 
when they are brought out of equilibrium by local pump¬ 
ing and extraction, thermally-activated collective effects 
inject a great amount of excitations in the system and 
drastically increase energy extraction from it. 

These results suggest to exploit the intrinsic non¬ 
locality of thermal noise as a resource to improve 
quantum energy transport, with direct applications to 
condensed-matter systems, to nanoscale energy manage¬ 


ment, and possibly being of great relevance also for light¬ 
harvesting 
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